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ABSTRACT 

We investigate the structural, kinematical, and spectrophotometric properties of 
“E+A” galaxies - those with strong Balmer absorption lines but no significant [Oil] 
emission - using numerical simulations combined with stellar population synthesis 
codes. We particularly focus on the two-dimensional (2D) distributions of line-of- 
sight velocity, velocity dispersion, colors, line index in E+A galaxies formed via the 
interaction and merging of two gas-rich spirals. Our numerical simulations demonstrate 
that E+A elliptical galaxies formed by major galaxy merging have positive radial color 
gradients and negative radial H(5 gradients by virtue of their central poststarburst 
populations. Furthermore, we show that the projected kinematical and spectroscopic 
properties of the simulated E+A galaxies can be remarkably different for different 
major merger models. For example, the simulated E+A ellipticals with kinematically 
decoupled cores clearly show regions of strong HJ absorption which are very flattened, 
with differences in rotation and velocity dispersion between the old and young stars. 
E+A ellipticals are highly likely to show more rapid rotation and a smaller central 
velocity dispersion in young stars than in old ones. E+A’s formed from the strong 
tidal interaction between gas-rich spirals have disky morphologies with thick disks 
and are highly likely to be morphologically classihed as barred SOs. We also provide 
specific predictions on the structural, kinematical, and spectrophotometric properties 
of young globular cluster systems in E+A’s. Based on these results, we discuss the 
importance of spatially resolved, integral field unit spectroscopy on large (8-lOm) 
ground-based telescopes in confirming the formation of kinematically distinct cores in 
elliptical galaxies produced via dissipative merging and determining the most probable 
physical mechanism(s) for E+A formation with disky and spheroidal morphologies. 

Key words: galaxies: starburst - galaxies:evolution - galaxies: elliptical and lentic¬ 
ular, cD - galaxies: kinematics and dynamics - galaxies: interaction 


1 INTRODUCTION 

Since Dressier & Gunn (1983) discovered the enigmatic 
”E+A” galaxies - i.e., those with strong Balmer line absorp¬ 
tion but no detectable emission - in distant clusters, their 
origin and the physical mechanisms responsible for their for¬ 
mation have long been discussed both observationally and 
theoretically. There have been two key issues which have 
driven much of this discussion: Firstly, what is responsi¬ 
ble for triggering and then abruptly truncating the recent 
starburst that most likely gave rise to the E+A spectral 
signature (Poggianti et al. 1999)? Secondly, what is the evo¬ 
lutionary link between E+A galaxies and the prevalent blue 
’Butcher-Oemler’ (Butcher & Oemler 1978) population? Al¬ 


though the nature of E+A galaxies in distant clusters has 
been extensively studied observationally (Couch & Sharpies 
1987, Barger et al. 1996, Couch et al. 1998, Dressier et al. 
1999, Galaz 2000, Tran et al. 2003), these two issues remain 
far from being resolved. 

Another important avenue of investigation has been to 
study the incidence of E+A galaxies in environments other 
than distant clusters, mainly at lower redshifts using large 
galaxy redshift surveys. Zabludoff et al. (1996, hereafter 
Z96) investigated 21 low-redshift (0.05 < 2 < 0.13) E+A’s 
identified in a sample of 11,113 galaxies observed as part of 
the Las Campanas Redshift Survey. They found that most 
of these E+A galaxies lay in the field rather than in clusters 
and therefore suggested that cluster environmental effects, 


© 1994 RAS 


2 K. Bekki, W. J. Couch, Y. Shioya, A. Vazdekis 


such as interactions with the cluster gravitational potential 
or intracluster medium, are not responsible for E+A forma¬ 
tion. More recently, Blake et al. (2004) conducted a simi¬ 
lar investigation using the order of magnitude larger 2dF 
Galaxy Redshift Survey (2dFGRS), and found the large ma¬ 
jority of the ^50-250 z ^ 0.1 E-l-A galaxies in their sample 
were either isolated or in galaxy groups, with only ~ 10 % 
being in rich clusters. At higher redshifts, Tran et al. (2004) 
investigated the fraction of E-l-A galaxies in a sample of ~ 
800 spectroscopically confirmed field galaxies in the range 
0.3 < 2 < 1.0 and found it to be ~ 2.7 ± 1.1%, which is sig¬ 
nificantly lower than that in galaxy clusters at comparable 
redshifts (11 ± 3%). 

Recent high resolution imaging, photometric and spec¬ 
troscopic studies of E-l-A galaxies undertaken with the 
Hubble Space Telescope (HST) and large ground-based 
telescopes have revealed a considerable diversity in the mor¬ 
phological, structural, and kinematical properties between 
E-|-A’s. Norton et al. (2001) investigated the internal kine¬ 
matics of the E-l-A galaxies in the Z96 sample and found 
that most of them were dynamically supported, with v/a 
and a ranging from 30kms“^ to 200kms“^ (here v is 
the rotational velocity and a is the velocity dispersion). A 
very recent HST-based morphological study of 5 of the Z96 
E-|-A’s by Yang et al. (2004), revealed that they were bulge- 
dominated systems (with and without a small underlying 
disk) and had radial luminosity profiles qualitatively similar 
to those of normal “power-law” early-type galaxies. They 
also discovered possible young and very bright (Mr < —13) 
globular cluster candidates around these E-l-A galaxies, with 
the number varying significantly from galaxy to galaxy. Tran 
et al. (2004) have also revealed the diversity in half-light ra¬ 
dius, total luminosity, morphological type, and internal ve¬ 
locity dispersion between E-fA’s in their samples. However, 
Blake et al. (2004) have cautioned that there maybe some 
heterogeneity and hence diversity in E-l-A galaxy samples 
due to the way they are selected; E-l-A samples whose strong 
Balmer line absorption has been identified only the basis of 
the equivalent width measured for the H(5 (and sometimes 
the H 7 ) line - as has generally been the case in distant clus¬ 
ter studies (Dressier et al. 1999) - will be contaminated by a 
population of disk-dominated, dusty star-forming galaxies, 
as well as containing bona fide E-|-A’s which are spheroid- 
dominated and have no ongoing star formation. 

This plethora of observational data on E-l-A galaxies, 
when taken in its entirety, raise many questions, the most 
significant being: (i) Gan model(s) of galaxy interactions 
and merging (e.g., minor vs major merging) explain self- 
consistently the observed dynamical, photometric, and spec¬ 
troscopic properties of E-|-A’s?. (ii) What mechanism drives 
the formation of disky E-|-A’s that dominates E-l-A popula¬ 
tions in clusters of galaxies? (iii) Are there any evolutionary 
links between the disky E-|-A’s in clusters and “passive spi¬ 
rals” with k-type spectra (Couch et al. 1998)? (iv) What is 
the origin of the positive and negative color gradients ob¬ 
served in E-fA’s? (v) Which can more self-consistently ex¬ 
plain the very red colors of E-fA’s discovered both in the 
field and in clusters - dust extinction or a truncated IMF? 
(vi) Are there any differences between the structural and 
kinematical properties of the old and young stellar popu¬ 
lations in E-l-A galaxies? (vii) Are the physical properties 
of bright young star clusters observed around E-|-A’s con¬ 


sistent with any formation scenario of globular cluster for¬ 
mation? Although previous theoretical studies have tried to 
understand possible star formation histories of E-l-A galaxies 
(Barbaro & Poggianti 1997; Poggianti et al. 1999; Shioya & 
Bekki 2000; Shioya et al. 2002; 2004), they had difficulty in 
addressing the above questions due to the limitations of the 
one-zone spectrophotometric models that they adopted. In 
order to properly tackle the above questions, the dynamical 
and spectrophotometric properties of E-l-A galaxies need to 
be investigated jointly in an explicitly self-consistent man¬ 
ner. Therefore, numerical simulations combined with stellar 
population synthesis codes - which enable us to predict not 
only structural and kinematical properties but also photo¬ 
metric and spectroscopic ones - are indispensable in solving 
the above important problems related to the formation and 
evolution of E-l-A galaxies. 

The purpose of this paper is to investigate the struc¬ 
tural, kinematical, chemical, photometric, and spectroscopic 
properties of E-fA galaxies in an explicitly self-consistent 
manner, thereby making considerable progress towards an¬ 
swering the above question. In particular, we investigate 
both radial gradients in the spectrophotometric properties 
as well as the projected (2D) distributions of Balmer ab¬ 
sorption, rotational velocity and velocity dispersion in E-l-A 
galaxies formed from galaxy interactions and merging: We 
will discuss other possible mechanisms of E-l-A galaxy for¬ 
mation in our forthcoming papers. These spatial proper¬ 
ties are the main focus of our study since current and 
planned spatially resolved spectroscopy with integral field 
units (IFUs) on 8 -lOm telescopes (e.g., GMOS on Gemini) 
is going to provide an unprecedented wealth of data on the 
two-dimensional distribution of spectroscopic properties for 
E-|-A’s (Pracey et al. 2004). Furthermore observational stud¬ 
ies based on the SAURON project with reasonably large 
field of views (e.g., Davies et al. 2001) now provide 2D spec¬ 
tra data of kinematics and stellar populations for the en¬ 
tire regions of galaxies. The simulated 2D distributions of 
the dynamical and spectroscopic properties will not only 
enable us to identify peculiar structural and kinematical be¬ 
haviour (e.g., kinematically decoupled cores) associated with 
the young stars in E-l-A galaxies, but also help us to cov¬ 
erage on the most reasonable E-l-A formation models. Our 
models can also predict the age-, metallicity- and spatial- 
distributions of young GCs associated with E-l-A galaxies, 
and so we can discuss in more depth the bright young can¬ 
didates recently observed by Yang et al. (2004). 

The plan of this paper is as follows: In the next sec¬ 
tion, we briefly describe our numerical models for E-l-A for¬ 
mation via galaxy interactions and merging. We summa¬ 
rize the numerical methods and techniques by which we 
can calculate the spectrophotometric properties of E-fA’s 
from N-body numerical data. In §3, we present the numer¬ 
ical results on the physical properties of E-|-A’s for differ¬ 
ent merger/interaction models. In §4, we use our results to 
address each of the seven questions mentioned above. We 
summarize our conclusions in §5. 

In presenting our work here, we stress that we are fo¬ 
cussing solely on those E-l-A galaxies formed as a result of 
galaxy-galaxy interactions and merging. This does not, how¬ 
ever, rule out other mechanisms (e.g., abrupt truncation of 
star formation via halo gas stripping and ram pressure strip¬ 
ping) being responsible for the formation of E-l-A galaxies. 
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Table 1. Model parameters 


model no. 

/b “ 

/g' 

m2 “ 

orbit 

comments 

Ml 

0.5 

0.1 

1.0 

PP 

major merger 

M2 

0.5 

0.2 

1.0 

PP 


M3 

0.5 

0.5 

1.0 

PP 


M4 

0.0 

0.1 

1.0 

IN 

no bulge 

M5 

0.5 

0.1 

1.0 

IN 

fiducial 

M6 

0.5 

0.2 

1.0 

IN 


M7 

0.5 

0.5 

1.0 

IN 


M8 

1.0 

0.1 

1.0 

IN 

more massive bulge 

M9 

0.0 

0.1 

1.0 

IN 

no bulge 

MIO 

0.5 

0.1 

1.0 

RR 


U1 

0.5 

0.1 

0.3 

IN 

unequal-mass merger 

U2 

0.5 

0.1 

0.1 

IN 


T1 

0.0 

0.1 

1.0 

PP 

tidal interaction 

T2 

0.0 

0.2 

1.0 

PP 


T3 

0.0 

0.5 

1.0 

PP 


T4 

0.0 

0.5 

1.0 

RR 


T6 

1.0 

0.1 

1.0 

PP 


T7 

0.0 

0.1 

0.1 

PP 


11 

0.0 

0.1 

- 

- 

isolated disk 

12 

0.5 

0.1 

- 

- 


13 

1.0 

0.1 

- 

- 



“ mass ratio of bulge to disk 
^ gas mass fraction 

mass ratio of two merging spirals 

and that different mechanisms can operate in different en¬ 
vironments (see, for example, Barbaro & Poggianti 1997; 
Bekki et al. 2001a,b; Bekki et al. 2002). 


2 THE MODEL 
2.1 Mergers 

Since the numerical methods and techniques we employ for 
modeling the chemodynamical and photometric evolution of 
galaxy mergers have already been described in detail else¬ 
where (Bekki & Shioya 1998, 1999), we give only a brief 
review here. 


2.1.1 Progenitor disk galaxies 

The progenitor disk galaxies that take part in a merger are 
taken to have a dark halo, a bulge, and a thin exponential 
disk. Their total mass and size are Ma and R^, respectively. 
From now on, all masses are measured in units of Ma and 
distances in units of i?a, unless otherwise specified. Velocity 
and time are measured in units of u = {GMa/Rdf^'^ and tayn 
= {Rfd/GMdf^'^, respectively, where G is the gravitational 
constant and assumed to be 1.0 in the present study. If we 
adopt Ma = 6.0 x 10^° M© and i?a = 17.5 kpc as fiducial 
values, then v = 1.21 x lO^kms”^ and tayn = 1.41 x 10® 

yr- 

We adopt the density distribution of the NFW halo 
(Navarro, Frenk & White 1996) suggested from CDM simu¬ 
lations: 

= (r/r,)(r+r/r.)2’ 


where r, po, and r^ are the spherical radius, the central den¬ 
sity of a dark halo, and the scale length of the halo, respec¬ 
tively. The value of rs (0.8) is chosen such that the rotation 
curve of the disk is reasonably consistent with observations. 
The bulge has a density profile with a shallow cusp (Hern- 
quist 1990): 

p{r) oc r~^{r + Obuige)”®, (2) 

where Ubuige is the scale length of the bulge and fixed at 0.04. 
The bulge mass and its compactness can control the bar for¬ 
mation in the disks and thus the strength of starbursts in 
mergers. The dark matter to disk mass ratio is fixed at 9 
whereas the bulge to disk ratio is assumed to be a free pa¬ 
rameter and represented by fh. The radial (R) and vertical 
(Z) density profiles of the disk are assumed to be propor¬ 
tional to exp(—7?/i?o) with scale length 7?o = 0.2 and to 
sec\A{Z/Zo) with scale length Zq — 0.04 in our units, respec¬ 
tively. In addition to the rotational velocity attributable to 
the gravitational field of the disk and halo components, the 
initial radial and azimuthal velocity dispersions are added to 
the disk component in accordance with the epicyclic theory, 
and with a Toomre parameter value of Q = 1.5 (Binney & 
Tremaine 1987) . The vertical velocity dispersion at a given 
radius is set to be 0.5 times as large as the radial veloc¬ 
ity dispersion at that point, as is consistent with the trend 
observed in the Milky Way (e.g., Wielen 1977). 

The disk is composed both of gas and stars, with the gas 
mass fraction (/g) being a free parameter and the gas disk 
represented by a collection of discrete gas clouds that follow 
the observed mass-size relationship (Larson 1981). All over¬ 
lapping pairs of gas clouds are made to collide with the same 
restitution coefficient of 0.5 (Hausman & Roberts 1984). The 
gas is converted into either field stars or globular clusters 
(GCs). Field star formation is modeled by converting the 
collisional gas particles into collisionless new stellar parti¬ 
cles according to the algorithm of star formation described 
below. We adopt the Schmidt law (Schmidt 1959) with ex¬ 
ponent 7 = 1.5 (1.0 < 7 < 2.0, Kennicutt 1998) as the con¬ 
trolling parameter of the rate of star formation. The amount 
of gas consumed by star formation for each gas particle in 
each time step is given by: 

Pg oc pf, (3) 

where pg is the gas density around each gas particle. The co¬ 
efficients in the law are taken from the work of Bekki (1998, 
1999) and the mean star formation rate in an isolated disk 
model for 1 Gyr evolution is ~ 1 M© for the adopted coef¬ 
ficient (thus consistent with observations). Globular cluster 
formation in the present model is discriminated from field 
star formation as follows. We use the cluster formation cri¬ 
teria derived by previous analytical works (e.g., Kumai et 
al. 1993), and hydrodynamical simulations with a variety of 
different parameters for cloud-cloud collisions on a 1-100 pc 
scale (Bekki et al. 2004) in order to model globular clus¬ 
ter formation. A gas particle is converted into a cluster if it 
collides with other high velocity gas (with the relative veloc¬ 
ities ranging from 30kms“^ to lOOkms”^) and having an 
impact parameter (normalized to the cloud radius) less than 
0.25. These stars formed from gas are called “new stars” (or 
“young stars”) whereas stars initially within a disk are called 
“old stars”, throughout this paper. 
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2.1.2 Chemical and spectrophotometric evolution 

Chemical enrichment through star formation during galaxy 
merging is assumed to proceed both locally and instanta¬ 
neously in the present study. We assign the metallicity of 
the original gas particles to the new stellar particles and in¬ 
crease the metallicity of the neighbouring gas particles. The 
total number of neighbouring gas particles is taken to be 
A^gas, given by the following equation for chemical enrich¬ 
ment: 

AMz = {^ii?metms-|-(1.0-7?met)(1.0-Zi)msJ/met}/Ai'gas.(4) 

Here, AMz represents the increase in metallicity for each 
gas particle, Zi the metallicity of the new stellar particle 
(or that of the original gas particle), 7?met the fraction of 
gas returned to the interstellar medium, nis the mass of the 
new star, and ymet the chemical yield. The values of i?met 
and j/met are set to 0.3 and 0.02, respectively. It is assumed 
here that the spectral energy distribution (SED) of a model 
galaxy is the sum of the SEDs of the individual stellar parti¬ 
cles. The SED of each stellar particle is assumed to be a sim¬ 
ple stellar population (SSP) that is a coeval and chemically 
homogeneous assembly of stars. Thus the monochromatic 
flux of a galaxy with age T, Fx{T), is described as 

F\{T) = Fssp^x{Zi,Ti) X nia , (5) 

star 

where Essp,A(^i, Vi) and iris are the monochromatic flux of a 
SSP of age Ti and metallicity Zi (where the suffix i identifies 
each stellar particle), and the mass of each stellar particle, 
respectively. The age of each SSP, Ti, is defined as n = T—ti, 
where ti is the time when a gas particle is converted into a 
stellar one. The metallicity of each SSP is exactly the same 
as that of the stellar particle, Zi, and the summation (Y.) 
in equation (5) is done for all stellar particles in a model 
galaxy. 

The adopted instantaneous recycling approximation 
means that the model does not consider the time delay be¬ 
tween star formation and the onset of supernovae explosion. 
This means that although chemical enrichment associated 
with type II SNe (that form ~ 10^ yr after star formation) 
can be modeled reasonably well, the one associated with 
SNela (that form ~ 10® yr after star formation) is not so 
correctly modeled owing to the large time delay between star 
formation and SNIa explosion. As a result of this, the present 
model can not correctly predict each chemical abundance 
(e.g.. Mg and Fe) in E-|-A’s: We have a plan to investigate 
chemical properties of E-|-A’s in our forthcoming papers by 
using more sophisticated models including SNIa effects. 

A stellar particle is assumed to be composed of stars 
whose age and metallicity are exactly the same as those of 
the stellar particle and the total mass of the stars is set to be 
the same as that of the stellar particle. Thus the monochro¬ 
matic flux of a SSP at a given wavelength is defined as 

/•Mu 

FsspAZi,Ti)= / 4>{M)fx{M,Ti,Zi)dM , (6) 

Jml 

where M is the mass of the star, and fx{M,Ti,Zi) is the 
monochromatic flux of a star with mass M, metallicity Zi 
and age n. 4>{M) is the initial mass function (IMF) of stars 
and Mu, Ml are the upper and lower mass limits of the IMF, 
respectively. In this paper, we use the Fssp,x{Zi,Ti) values 


calculated by Vazdekis et al. (1996) with an IMF slope of 1.3 
(i.e., the Salpeter IMF). We also adopt their two different 
models with Ml = 0.1 M© and 1.0 Mq. 

In order to calculate the spectrophotometric evolution 
of merger remnants with poststarburst components, we need 
to assume that (l)the initial stellar disk and bulge in a 
merger progenitor spiral have the same age of Taiak, and 
(2) the disk has a metallicity gradient consistent with obser¬ 
vations. Considering the recent observations of the Galactic 
metallicity gradient (Friel 1995), we allocate metallicity to 
each disc star according to its initial position: At r = R, 
where r{R) is the projected distance (in units of kpc) from 
the center of the disc, the metallicity of the star is given by: 

[m/H]^^j^ = [m/H]^^^^P -b Od x R. (7) 

We adopt a plausible value of —0.091 for the slope aa from 
Friel (1995), in which the Galactic stellar metallicity gra¬ 
dient is estimated from Galactic open clusters. The central 
value of [m/H]j is chosen such that [m/H]j,^j^ is 0.02 at 
the radius of 8.5 kpc (the solar radius for the Galactic disk). 
The merger models with the adopted metallicity gradients 
above are demonstrated to explain consistently the observed 
metallicity distribution function of the stellar halo in NGG 
5128 (Bekki et al. 2003). We investigate the models with 
Tdisk = 5, 7, 10 Gyr, because our results can be compared 
with the relatively nearby (2 ~ 0.1 — 0.3) E-|-A’s. We show 
only the results with Tdisk = 7 Gyr in the present study, 
because the results do not strongly depend on Tdisk for the 
above range of Tdisk. Several authors have suggested that 
starbursts produce a top-heavy IMF, because the formation 
of lower-mass stars is suppressed during such an event (e.g., 
Larson 1998). Gonsidering this possible difference in the IMF 
between starburst phases and periods of normal star forma¬ 
tion, we adopt Ml =0.1 M© for old stars and Ml = 1.0 
Mq for new stars formed during starbursts in deriving Tssp 
for each stellar component. 

The main reason for our adopting a ’top-heavy’ IMF for 
the young starburst component is as follows: Shioya et al. 
(2004) have demonstrated that poststarburst models with 
a normal IMF (with Ml = 0.1 Mq) are unable to simul¬ 
taneously explain the red colors and strong H5 absorption 
that are seen for a subset of the E-l-A galaxies observed by 
Gouch & Sharpies (1987), Balogh et al. (1997), and Galdwell 
et al. (1999). We need either dust extinction (e.g., Shioya et 
al. 2004) or a top-heavy IMF (Balogh et al. 1997; Ghar- 
lot et al. 1993) to explain these red-Hd-strong E-|-A’s where 
starbursts have ended: A top-heavy IMF is a possibility for 
explaining the origin of these red E-l-A galaxies. Therefore 
we have chosen a top-heavy IMF for the young stars in the 
present study. If the normal IMF is adopted for young stars, 
the poststarburst signature becomes rather weak in E-l-A 
phases. For example, EW(H5) in the standard model is 7.4 
A for Ml = 1.0 M© and 1.1 A for Ml = 0.1 M© at the 
same time T. 

The present model does not consider the effects of dust 
extinction on the SEDs of poststarburst galaxies, nor the 
effects of gaseous emission on the Balmer absorption lines 
(It should be stressed here that spectral indices are not af¬ 
fected by dust emission and Balmer lines might be affected 
by nebular emission). Since most of the gas in an interact¬ 
ing/merging pair should be consumed in the associated star- 
burst and hence by the time it develops an E-l-A spectral 
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signature, the omission of dust extinction effects is justified. 
We nonetheless discuss the possible important effects of dust 
in E+A spectra in §4. In order to incorporate gaseous emis¬ 
sion into our model SEDs, we need to combine our N-body 
simulation code with both the adopted SSP code (Vazdekis 
et al. 1996) and with one that can provide emission lines 
luminosities (e.g., PEGASE by Eioc & Rocca-Volmerange 
1997). To combine these three different numerical codes is a 
formidable task, particularly into our new chemodynamical 
model which already includes a variable IMF and globular 
cluster formation. Hence we will leave a full treatment of 
gaseous emission for a future paper. We therefore adopt the 
assumption that the dilution of the Balmer absorption fea¬ 
tures by emission from interstellar gas ionized by residual 
star formation in the post-starburst phase does not occur 
at all. Star formation is therefore assumed to be completely 
halted ~ 1.0 Gyr after the star formation associated with 
the starburst reaches its maximum rate. We therefore do 
not intend to discuss the origin of a small amount of possi¬ 
bly residual star formation observed in a few E+A galaxies 
(Miller & Owen 2001). AGN feedback effects from massive 
black holes on spectrophotometric properties of E+A’s can 
be also important (Shioya & Bekki 2005). 

2.1.3 Orbital configurations 

In all of the simulations of merging pairs, the orbit of the two 
disks is set to be initially in the xy plane and the distance 
between the center of mass of the two disks is assumed to be 
10 in our units (corresponding to 175 kpc). The pericenter 
distance and the eccentricity are set to be 1.0 (17.5 kpc) and 
1.0 (i.e., parabolic), respectively, for most of the models. The 
spin of each galaxy in a merger is specified by two angles 6i 
and (pi, where suffix i is used to identify each galaxy. 6i is 
the angle between the 2 axis and the vector of the angular 
momentum of a disk, fii is the azimuthal angle measured 
from the x axis to the projection of the angular momentum 
vector of a disk onto the xy plane. We specifically present 
the results of the following three models with different disk 
inclinations with respect to the orbital plane: A prograde- 
prograde model represented by “PP” with 9\ = 0, 02 = 30, 
(()i = 0; a retrograde-retrograde model (“RR”) with 9\ = 
180, 02 = 210 , (j)\ = 0 , and <j )2 ~ 0 ; and a highly inclined 
model (“IN”) with 0i = 60, 02 = 6Q, 4>\ =90, and (f >2 = 
0. The time taken for the progenitor disks to completely 
merge and reach dynamical equilibrium is less than 16.0 in 
our units (~ 2.2 Gyr) for most of our major merger models. 
We also present the results of unequal-mass merger mod¬ 
els with the mass ratios of the two merging spirals (repre¬ 
sented by m 2 ) equal to 0.1 and 0.3. Table 1 summarises the 
model parameters for the merger models, with major merger 
models and unequal-mass models labeled as “M” and “U”, 
respectively. 


2.2 Tidal interaction 

The initial disk models and the numerical methods used for 
computing chemodynamical and spectrophotometric evolu¬ 
tion are exactly the same for the tidal interaction models 
as they are for the merger models. In the tidal interaction 
models, the two disk galaxies do not merge with each other 
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Time (Gyr) 

Figure 1. Star formation histories of major galaxy mergers be¬ 
tween gas-rich spirals for the model M2 with /g = 0.2 (thin solid) 
and the model M5 with /g = 0.5 (thick solid) for 2.4 Gyr. These 
two models are examples of merger-driven massive starbursts with 
a subsequent rapid decline in star formation rate. 

this possible difference of IMF between starburst phases and nor¬ 
mal ones in galaxy evolution, 
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Figure 2. Mass distribution of stars of the simulated E+A’s pro¬ 
jected onto the x-z plane (corresponding to the orbital plane) for 
four representative models: Ml (upper left), M5 (fiducial, upper 
right), U1 (lower left), and T1 (lower right). EW(H<5) in each 
model is shown in the upper part of each frame. Note that all 
models show large EW(H<5) (> 6 A). 

by virtue of the large pericenter distance (2 in our units, 
corresponding to 35 kpc) that is adopted. Therefore, only 
one of the two interacting galaxies needs to be modeled as 
a fully self-gravitating N-body system, while the other is 
modeled as a point-mass particle. We mainly present the re¬ 
sults of the tidal interaction models with m 2 = 1.0: Strong 
starbursts do not happen in interaction models with smaller 
m 2 (~ 0 . 1 ); hence these models do not have poststarburst 
populations with strong E+A spectra. Moderately strong 
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starbursts, with a maximum star formation rate of 20 — 30 
Mq, occur in tidal interactions with m2 = 1.0; consequently 
most of the gas within interacting disks is consumed within 
~ 1 Gyr. 

Table 1 summarises the model parameters for the tidal 
interaction models, all of which are labeled “T”. For com¬ 
parison, we also investigated isolated disk models (without 
tidal interaction and merging) and these models are labeled 
as “I” in Table 1. The results of these isolated models are 
used when we investigate how strongly star formation rates 
are increased in interacting/merging galaxies compared with 
an isolated disk with no external perturbation. Star forma¬ 
tion rates in these isolated disk models are at most a few 
Mq yr“^ whereas those in interacting/merging models are 
an order of 10 — 100 Mq yr“^. The results of the isolated 
models are only very briefly described in the present study, 
simply because they do not show E-l-A spectra. 

All the calculations related to the above chemodynam- 
ical evolution have been carried out on the GRAPE board 
(Sugimoto et al. 1990) at the Astronomical Data Analysis 
Genter (ADAG) at the National Astronomical Observatory 
of Japan. The gravitational softening parameter was fixed 
at 0.025 in our units (0.44 kpc). The time integration of the 
equation of motion was performed by using the 2nd-order 
leap-frog method. The initial total particle number in each 
simulation was 110,000 for a ‘merger’ model and 55,000 for 
the ‘interaction’ and ’isolated’ models. 

2.3 Main points of analysis 

Star formation histories and morphological evolution of 
interacting and merging galaxies with strong starbursts 
have already been described in previous numerical studies 
(e.g., Noguchi & Ishibashi 1988; Mihos & Hernquist 1996; 
Bekki 1998). Furthermore, spectral type evolution (e.g., 
from “e(a)” to “a+k”; Dressier et al. 1999) in starburst¬ 
ing galaxy mergers has been extensively discussed by Bekki 
et al. (2001a) based on spectro-dynamical simulations. In 
this paper, therefore, we particularly focus on (1) the spatial 
distribution of the optical and near-infrared colors and the 
equivalent widths of the Balmer absorption lines in E-|-A’s, 
and (2) the difference in the kinematical properties of the 
old and young stellar components in E-l-A galaxies. 

As far as (2) is concerned, Norton et al. (2001) found 
possible differences in rotational velocity and velocity dis¬ 
persion between the old and young components in nearby 
E-|-A’s. Our investigation of (2) therefore involves checking 
whether the simulated kinematics of E-l-A galaxies formed 
by galaxy interactions and merging are consistent with such 
observations. The kinematical properties of the young stel¬ 
lar populations in E-|-A’s are derived from the properties 
of the strong Balmer absorption lines, which are the key 
observables in this context. To compare these observables 
with the present numerical results in a consistent way, we 
show the kinematics of young H((5) strong stars that have 
EW(H(5) larger than 2A. We attempt here to completely re¬ 
move the contribution from new stars that are formed in 
the very early phase of merging (and thus have relatively 
old ages) and consequently do not contribute to the strong 
Balmer line absorption seen in the E-l-A phase. 

We also investigate the physical properties of the simu¬ 
lated GG’s in E-|-A’s, because we consider that the predicted 


properties of young GC’s in E-|-A’s are useful in discussing 
whether the adopted cloud-cloud collision model of GG for¬ 
mation is consistent with the observations of young blue GG 
candidates in Ej-A’s (e.g., Yang et al. 2004), and whether 
there are evolutionary links between blue GG candidates 
in E-|-A’s and the red metal-rich GG’s around passive el¬ 
lipticals. We concentrate in particular on the 2D dynamical 
and spectroscopic distributions and the GG properties of our 
model merger/interaction remnant galaxies when they have 
an HS equivalent width, EW(H5), larger than hA, consistent 
with stringent criterion adopted by Z96 in their E-l-A selec¬ 
tion. We here stress that when the simulated model shows 
strong EW(HJ), other Balmer absorption lines also show 
large values: The simulated models with EW(HJ) larger than 
6A are consistent with Z96’s selection criterion. 

The star formation histories, which determine the fi¬ 
nal spectrophotometric properties of the remnant interact¬ 
ing/merging galaxies, differ quite significantly in our simu¬ 
lations, being dependent on the gas mass fraction (/g), the 
bulge-mass fraction (/b), and the mass ratio of the merg¬ 
ing two disks (m2). Figure 1 shows an example of this de¬ 
pendence. The structural and kinematical properties of the 
simulated E-|-A’s are also diverse, depending mainly on m2 
and the orbital configurations. We therefore describe the re¬ 
sults of our most representative models which show the typi¬ 
cal and/or most interesting behavior in E-l-A formation and 
evolution. We first describe the physical properties of the 
simulated E-l-A in the fiducial merger model M5, and then 
discuss their parameter dependences based on the results 
of other representative merger models. The morphological 
properties seen from the orbital plane of the galaxy interac¬ 
tion/merger (e.g., the x-z plane) are shown in Figure 2 for 
the representative models, Ml, M5, Ul, and Tl. We mainly 
describe the results of these representative models. 

We will explain separately the results of the merger 
models and the tidal interaction models, in order to con¬ 
trast the differences in properties between the two. In the 
following sections, the time T represents the time that has 
elapsed since the simulation starts. 

3 RESULTS 

3.1 The fiducial merger model 

3.1.1 2D distributions of speetrophotometric properties 

Figure 3 summarises the global morphological properties of 
the stellar populations with different EW(H(5) for the E-l-A 
formed in the fiducial model at T = 2.8 Gyr. In this major 
merger model with a highly inclined orbital configuration 
and a large pericenter distance (17.5 kpc), the old stars with 
EW(HJ) < 0 A appear to have the morphology on an el¬ 
liptical galaxy. Furthermore, there is no sign of any tidal 
features, because the tidal tail formed during the merger 
disappeared by the time the remnant shows strong Balmer 
line absorption. The major axis of the E-l-A is nearly coin¬ 
cident with the x-axis for the old stars. 

The young stars that are formed during the merger 
(with 0 <EW(H(5) < 6A), have more compact spatial dis¬ 
tributions than the old stars, although their major axes are 
nearly aligned with those of the old stars. Young stars with 
EW(H(5) > 6A have the most compact spatial distribution. 
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Figure 3. Mass distributions projected onto the x-z plane in 
the fiducial model at T = 2.8 Gyr for stellar populations with 
different EW(H(5): For EW(H(5) < 0 A (upper left), 0 < EW(H5) 

< 2 A (upper right), 2 < EW(H<5) < 6 A (lower left), and 6 A 

< EW(H(5) (lower left). 
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Figure 5. The two-dimensional (2D) distributions of EW(H5) 
(left) and EW(H/3) [right) projected onto the x-z plane in the 
fiducial model at T = 2.8 Gyr. The abscissa and the ordinate 
represent the x-axis and the 2 :-axis, respectively. Here we divide 
the 3 kpc X 3 kpc central region (shown in Fig. 4) into 20 x 20 
grid points and thereby estimated the SED of each grid point. For 
clarity, the grid points with EW(H(5) (and EW(H/3)) less than OA 
are shown in the darkest color. Note that the 2D distributions 
show the strong Balmer line absorption to be elongated along the 
2 -axis. 
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Figure 4. The same as Fig. 3 but for the central 3 kpc. 


and are confined to the central few kiloparsecs. This reflects 
the fact that these stellar populations were formed in the 
last major starburst during the merger. 

Figure 4 shows that there is a clear difference in 
the morphological properties between the old stars (with 
EW(H(5) < OA) and the youngest stars (with EW(H(5) > 
gA) within the central 3 kpc of the E+A formed in the fidu¬ 
cial model. In the central 3 kpc, the major axis of the old 
stars’ mass distribution is nearly aligned with the x-axis 
and thus with that of the global mass distribution (shown 
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Figure 6. The same as Fig. 5 but for the optical color V — I (left) 
and the near-infrared color I — K [right). 
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Figure 7. Radial profiles of F — 7 [lef£) and EW(H5) [right) 
derived from V — I and EW(H(5) of the 400 grid points shown in 
Figs. 5 and 6. 
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Figure 8. The same as Fig. 5 but for T = 4.5 Gyr, when the 
model shows a passive (“k”-type) spectrum. 


in Fig. 1), whereas the major axis of the young stars is 
nearly parallel to the 2-axis (and thus perpendicular to that 
of the mass distribution of old stars). This misalignment in 
the mass distributions of the E-l-A results from the dissipa¬ 
tive formation of a thick stellar disk composed of stars with 
strong Balmer absorption lines in the E-l-A core. The more 
diffuse distribution of young stars (with 0 <EW(H5) < SA), 
is due to the earlier formation of these stars within the 
merger progenitor disks (rather than in the very center of 
the merger at the final phase of merging). 

Figure 5 clearly shows that the 2D distributions of 
EW(H5) and EW(H/3) absorption have very flattened shapes 
along the 2-axis, and the direction of elongation is coincident 
with the major axis of the mass distribution of young stars 
with EW(H(5) > 6A(shown in Fig. 4). It also shows that both 
EW(H5) and EW(H/3) are larger in the inner regions of the 
E-l-A, and therefore have a negative radial gradients. These 
results are all due principally to the centralized starbursts 
during dissipative merging in this model. The flattened 2D 
distribution can also be seen in the optical {V — I) and near- 
infrared {I — K) color distributions shown in Figure 6, which 
confirms that the flattened shape is due to the flattened dis¬ 
tribution of young stars in the core of the E-|-A’s. 

Figure 7 shows that the radial gradient m V — I color 
is “positive” in the sense that the inner color is bluer; this 
is also true for I — K. This figure also shows graphically, 
the negative radial gradients seen in EW(H(5) and EW(H/I). 
We find that the mean logarithmic gradient m V — I, i.e., 
A(I7 —/)/A(logi?) is O.ldmag per dex in radius, R, whereas 
that in EW(H(5) (A775/A(logi?)) is —2.6 A per dex in ra¬ 
dius R. These positive color gradients and negative EW(Hd) 
gradients in our simulated E-|-A’s are consistent with what is 
observed in some E-|-A’s (e.g., Yang et al. 2004). We discuss 
such comparisons further in §4. 

In the E-l-A simulated in the fiducial model, the spec- 
trophotometric properties - as delineated by the 2D distri¬ 
butions - rapidly evolve with time, due to the fading of the 
young stars. Figure 8 shows that the elongated distribution 
in Balmer line absorption seen in the strong E-l-A phase 
(T = 2.8 Gyr) becomes nearly invisible at T = 4.5 Gyr (i.e., 
1.7Gyr after the strong E-l-A phase). This disappearance 



Figure 9. Distributions in the number of grid points with a given 
color {left) and EW(H5) (right) for two different epochs, T = 
2.8 Gyr (thick solid line) and T = 4.0 Gyr (thin solid line), in the 
fiducial model. The color and EW(H5) values are those at each 
grid point in the 2D distributions (as shown in Figs. 5 and 6) and 
the normalized number of grid points is shown for convenience. 
The spectral type of the simulated E-l-A evolves from “a+k” to 
“k-|-a” (to ”k”) during these two epochs. Note that the dispersion 
of the distribution becomes very small in the k+a phase (i.e., 
T = 4.0 Gyr) both for F - / and EW(H5). 
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Figure 10. The 2D distribution of line-of-sight velocities viewed 
from the y-axis (i.e., line-of-sight velocity field projected onto the 
x-z plane) for old stars (left) and young stars (right) in the fiducial 
model at T = 2.8 Gyr. The abscissa and the ordinate represent the 
a;-axis and the 2 -axis, respectively. This figure enables us to con¬ 
firm whether the simulated E+A has global rotation and which 
direction the ZVC extends. The ZVG extends from the lower left 
corner to the upper right corner in this figure for old stars. Each 
frame measures 6kpc so that this figure can be compared with 
Figs. 4, 5, and 6. 


of the flattened shapes is also seen in the 2D V — I and 
I — K distributions. Figure 9 describes how the 2D V — I 
and EW(Hd) distributions evolve from the E-l-A phase to the 
passive (”k” type) phase. It is clear from this figure that the 
dispersion in the distribution becomes significantly smaller 
as the E-l-A elliptical evolves into a passive system both for 
F - / and EW(H5). The F - / (EW(H5)) dispersion de¬ 
creases from 0.12 mag (2.2 A) to 0.04 mag (0.8 A) during the 
transition phase. These results suggest that the dispersion in 
the 2D spectrophotometric properties can be regarded as an 
evolutionary “time arrow” in going from the E-l-A phase to 
passive phase. We thus suggest that elliptical galaxies with 
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Figure 11. Radial profiles of line-of-sight velocity along the 2 - 
axis derived from grid points shown in Fig. 10 for old stars (left) 
and young stars [right). Note that the radial gradient is steeper for 
the young populations than for the old ones, though the dispersion 
in the outer part is large for the young ones. 
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a of young stars 
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Figure 12. The same as Fig. 10 but for the velocity dispersion. 
Note that the young stars show smaller dispersion in the inner 
regions of the E+A. 


weaker H(5 have smaller dispersion in the 2D distributions 
of colors and Balmer absorption lines. 


3.1.2 Kinematics 

Figure 10 displays the 2D distribution of line-of-sight- 
velocities for the old and young stars for the simulated E+A 
of the fiducial model at T = 2.8 Gyr. It can bee seen that 
for the old stars, the zero velocity curve (ZVC) - which is 
defined as a line connecting the points where the line-of- 
sight velocity is zero - extends from the lower left corner 
to the upper right corner and thus is not aligned with the 
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Figure 13. The same as Fig. 11 but for the velocity dispersion. 


major axis of the mass distribution of old stars shown in 
Figs. 3 and 4. The direction of the ZVC of young stars with 
EW(H5) > 2 A, is broadly consistent with that of the ZVC 
of old stars, though there is a clear difference in the veloc¬ 
ity field along the a-axis (i.e., the vertical axis) between the 
two components for \x\ < 2kpc. The misalignment seen be¬ 
tween the ZVC and the major axis of the mass distribution 
in E+A’s, both for old and young stars, clearly indicates the 
minor-axis of rotation of the E+A. 

Eigure 11 shows the radial rotation profile along the 
minor axis (i.e., 2-axis) derived from the data points shown 
in Figure 10. Figure 10 shows that the rotation profile along 
the minor axis is significantly different between the old stars 
and the young stars. By applying a least-squares fit method 
to the grid points with |®| < 1.5 kpc in Figure 10, we can 
quantify the difference in the minor-axis rotation between 
the two. The derived rotational velocity, Kot, expressed as 
a function of the 2 coordinate (in units of kpc), is given by: 

Kot «-18.5 X 2-2.5 (8) 

for the old stars, and 

Kot «-27.0 X 2-1.2, (9) 

for the young stars, where the velocities are in units of 
kms“^. In the same way, the major axis rotation profile 
(along the ®-axis) can be estimated. The derived Kot, ex¬ 
pressed as a function of the x coordinate is: 

Kot « 4.4 X a; - 1.9 (10) 

for the old stars, and 

Kot « 9.3 X a; - 0.2 (11) 

for the young stars. 

These results clearly indicate that the young stars are 
more strongly supported by rotation than the old stars in 
the model E+A. The origin of the steeper radial gradient of 
Kot along the minor axis is closely associated with the for¬ 
mation of the kinematically distinct core (KDC) composed 
mostly of young stars with strong H5 absorption in the cen¬ 
tral few kpc of the E+A. Owing to the efficient gaseous dis¬ 
sipation during the very late phase of galaxy merging, the 
central core can have a significant amount of rotation with 
the angular momentum axis different from that of the ma¬ 
jor old component. As a result of this formation of a KDC, 
the minor axis rotation can be more remarkable in young 
stars than in old ones in the 2D velocity distribution. Thus 
Figures 5 and 10 suggest that future observations can con¬ 
firm the dissipative formation of KDCs in galaxy mergers by 
investigating the 2D distribution of line-of-sight velocities of 
young stellar populations in E-hA’s. 

Figure 12 shows the 2D distribution of line-of-sight ve¬ 
locity dispersion (ct) both for old and young stars in the 
simulated E+A. Although the dispersion field appears to 
be somewhat irregular, in particular, for the young stars, it 
shows that the central dispersion is higher for the old stars 
than the young stars and the radial gradient of the disper¬ 
sion is “positive” for young stars in the sense that the inner 
velocity dispersion is lower in the inner regions than in the 
outer regions for the central few kpc of the E+A. These two 
results do not depend on parameters of the present merger 
models, so that they can be considered to be generic char¬ 
acteristics of the simulated E+A’s. 
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Figure 14. The distributions of young GC’s with EW(H<5) > 6 A 
(upper two panels), and with EW(H5) < bA (lower two panels), 
projected onto the x-y plane (left two panels) and onto the x-z 
plane (right two panels) in the fiducial model with T = 2.8 Gyr. 


Figure 13 shows that the radial gradient of velocity dis¬ 
persion is steeper for the young stars the older ones. If we 
assume a oc R, where R is the distance of each grid point in 
Figure 12 from the center of the E-l-A, then we can roughly 
quantify the difference in the radial gradient between the 
two components based on the least-squares fit method. The 
variation of the derived velocity dispersion, a, with radius, 
R, for R < 3 kpc is described by: 

a «-11.5 X i?-h 210.0 (12) 

for the old stars, and 

a « 8.3 X i?-h 130.0 (13) 

for the young stars. The reason for the young stars appar¬ 
ently having such a small a slope is that the dispersion in a 
is large both at larger radii (2 — 3 kpc) and at the centrer 
(See Fig. 13): The least-squares fit to data points with the 
projected distance (R) less than 3 kpc give a small a slope 
owing to this large dispersion. It should be stressed here that 
if the cr slope in young stars is estimated for data points with 
R less than 1.5 kpc, the slope becomes significantly steeper 
(cr oc 35 X i?). We discuss the consistency of these results 
with recent observations by Norton et al. (2001) later in §4. 

It should be stressed here that the positive radial gradi¬ 
ent in velocity dispersion we see for the young stars, is con¬ 
sistent with the radial velocity dispersion profiles observed 
in bright nucleated dwarf elliptical galaxies (dE,Ns) in the 
Virgo cluster (Geha et al. 2002). This consistency implies 
that dissipative processes associated with nuclear starbursts 
can be responsible for the formation of the lower velocity dis¬ 
persion nuclei in dE,N’s. Given the fact that the velocity dis¬ 
persion of the young stars in the E-l-A galaxies, EA5, EA7, 
and EA15, observed by Norton et al. (2001) also decreases 
inwardly, there could be an evolutionary link between such 
low luminosity E-|-A’s and dE,Ns. 
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Figure 15. The distributions of young GC’s on the (V — /)- 
My plane (upper left), the (I — K)-My^ plane (upper right), the 
(B — V)-(V — I) plane (lower left), and the (B — U)-EW(H5) plane 
(lower right), for the fiducial model with T = 2.8 Gyr. 


3.1.3 GC properties 

Figure 14 shows the spatial distributions of young GG’s 
formed from high-speed cloud-cloud collisions with small im¬ 
pact parameters in the fiducial model. The mass fraction of 
GC’s to field stars in this model is 7%, that is a factor of ~8 
larger than that in the isolated model II (1%). This suggests 
that GC formation is more strongly enhanced than field star 
formation in dissipative major merging. Nearly all (100 %) of 
young GC’s with H5 > 6 A are within the central 3 kpc (i.e., 
the projected distance of R) less than 3 kpc) of the E-|-A 
whereas about 70% of young GCs with H5 < 6 A are within 
the central 3 kpc, which indicates an age-dependent spatial 
distribution for the young GC’s in the E-l-A (i.e., a more 
compact distribution for the younger GC’s). This result im¬ 
plies that the vast majority of Hd strong GC’s will be very 
difficult to identify, observationally, as GC’s associated with 
the E-l-A, due to their close proximity to their host’s center. 
Only the outer young GC’s with relatively small EW(H(5) 
will be detected, appearing as relatively bright blue com¬ 
pact sources. Most of these young GC’s have a metallicity 
that is higher than solar, and hence will evolve into red, 
metal-rich GC’s a few Gyr later. 

Figure 15 shows that young GC’s in the E-l-A have rel¬ 
atively bright luminosities with My >-10 mag and Mk > 
-12 mag. This is because of the young age of the GC’s and 
the relatively high value (1.0 M©) that we have adopted for 
the lower-mass cut off of the IMF (mu). It should be noted 
here that if we adopt the lower value of 0.1 M© for mu, the 
above absolute magnitudes of the GC’s become significantly 
(several magnitude) fainter. Since young GCS are all as¬ 
sumed to have the same mass in the present study, the older 
(fainter) GC’s will have redder colors, as shown in Figure 15. 
The dispersion in the optical colors, V — I, oi the GC’s is 
large, not because the GC’s have diverse metallicities but 
because there are age differences between them. The GC’s 
are actually distributed with a small dispersion along a line 
in the (B — V)-(V — 1) plane, due to their narrow range in 
metallicity. 
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Figure 16. The same as Fig. 10 but for the model Ml at T = 
2.8 Gyr. 
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Figure 17. The same as Fig. 10 but for the model U1 at T = 
2.8 Gyr. 


The locations of the GC’s on the {B - R)-EW{HS) 
plane in Figure 15 clearly demonstrates that they are on the 
“poststarburst evolutionary path” first identihed by Couch 
& Sharpies (1987), where they evolve from blue H(5 strong 
objects to red H(5 weak objects. Hence blue strong and 
red H(5 weak GC’s appear to coexist in E+A’s formed via 
major merging with strong starbursts. This results also im¬ 
plies that E+A’s formed via other physical mechanisms, 
such as truncation of star formation without a starburst, do 
not show such coexistence in their GC systems. Red GC’s 
in passive ellipticals can be either GC’s that were initially 
associated with the merger progenitor spiral or those that 
were formed during dissipative merging. Bright, blue GC’s 
with strong H(5 absorption lines in E+A’s are unambiguously 
young GC’s formed during dissipative merging. We therefore 
suggest that if the derived age-dependent spatial distribu¬ 
tion and spectrophotometric properties of young GC’s in the 
simulated E+A’s can be compared with future observations, 
such theories on the formation of GC’s during mergers can 
be tested in a more stringent way. 


3.2 Parameter dependences of merger models 

3.2.1 Generie results 

Before discussing the dependencies our merger model results 
have on the various physical parameters, it is useful to sum¬ 
marise the more generic results of our study. These are as 
follows: 

(i) E+A’s formed by galaxy merging exhibit positive 
color gradients and negative EW(H5) gradients, due to the 
poststarburst populations that reside within a few kpc of 
their centers. By mass, the fraction of all young stars with 
EW(H(5) > 2 A among all new stars can be as high as ^60% 
in the central 3 kpc of these E+A’s. The mass fraction of 
young stars with EW(H(5) > 2 A among all stars (i.e., old + 
new stars) can become as high as 0.2 in the central 3 kpc. 
This fraction depends on the gas mass ratio of /g in such 
a way that the fraction is larger for larger /g (e.g., 0.16 for 
fg = 0.1 and 0.20 for fg = 0.5). The fractional light from 
these poststarburst populations in the standard model for 
B—band is ~ 60 % for Ml =0.1 M© and ^ 100 % for 1.0 M© 
at the strong poststarburst epoch: The SEDs of E+A’s can 
be largely determined by young stars. 

(ii) These young stars (with EW(H(5) > 2 A) show more 
rapid rotation and a smaller central velocity dispersion than 
the old stars in the central 3 kpc of E+A’s. The radial gra¬ 
dient in the velocity dispersion of these stars is positive in 
the sense that the velocity dispersion is smaller in the inner 
regions of an E+A. 

(iii) E+A’s formed in this way also have young, H5- 
strong, metal-rich (more than the solar metallicity) and blue 
GG’s, most of which are located in their central regions. The 
GG formation efficiency increases dramatically (by a factor 
of 10) during dissipative merging, so that E+A’s end up 
with a larger GG specific frequency (i.e., number of GGs per 
unit luminosity) compared to their progenitor disk galaxies. 
Number fraction of GG particles to all new stellar ones is ~ 
1 % for the isolated disk models, II, 12 and 13 and ~ 7 % 
for the major merger model Ml. 

3.2.2 Orbital configurations 

Another key feature of our models is that they show the 
kinematical properties of the old and young stars in E+A’s 
depend strongly on the orbital configurations of the major 
merger event responsible for their formation. This depen¬ 
dency can be summarised as follows: 

(i) Some E+A ellipticals show a clear sign of major-axis 
rotation, yet no sign of minor-axis rotation in the 2D dis¬ 
tribution of line-of-sight velocities for both old and young 
stars. Figure 16 shows an example of this for the major 
merger model Ml in which an E+A elliptical is formed via 
a nearly “prograde-prograde” merger, with the orbital spin 
axis nearly parallel to the intrinsic spin axes of the two 
merger progenitor disks. The rotational velocity, I+ot (km 
s“^), expressed as a function of the x coordinate (kpc), that 
is derived for the 2D velocity distribution of the young stars 
in this model for \x\ < 3kpc is: 

Kot « 20.1 X a: +13.9 (14) 

for old stars and 

Kot « 21.7 X a: + 13.0. (15) 


© 1994 RAS, MNRAS 000, 1-?? 




12 K. Bekki, W. J. Couch, Y. Shioya, A. Vazdekis 


o of old stars a of young stors 



-2 0 2 -2 0 2 

kpc kpc 


0 50 100 150 200 0 50 100 150 200 

km/s km/s 


Figure 18. The same as Fig. 11 but for the model Ul. 
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Figure 19. The same as Fig. 13 but for the model Ul. 

The major-axis rotation is much less remarkable for the 
model MIO in which an E-l-A is formed via a “retrograde- 
retrograde” merger. 

(ii) Quite flattened EW(H(5), EW(H/3), V-I, and I-K 
distributions are formed in the central few kpc of the E-|-A 
if the major merger responsible has a highly inclined orbital 
configuration and at least one disk orbiting in a retrograde 
manner (i.e., the intrinsic spin axis is anti-parallel to its 
orbit with respect to the mass center of the merger). The 
origin of the flattened distributions is due essentially to the 
formation of KDCs composed of young, HM-strong stars. Due 
to the presence of young KDCs, E-l-A ellipticals can show 
a significant difference in the 2D line-of-sight velocity and 
velocity dispersion between their old and young stars. 

3.2.3 Mass ratio, m 2 

The dependences of E-l-A properties on the merging galaxy 
mass ratio, m2, can be described as follows: 

(i) E-|-A’s formed by unequal-mass mergers are highly 
likely to show rapid rotation along their major-axis, irre¬ 
spective of the orbital configuration of the merger, if they 
are viewed edge-on. Figure 17 shows the 2D distributions 
of the line-of-sight velocities for old and young stars in the 
model Ul, where more rapid rotation is seen than in the 
Ml model shown in Figure 16. This is due to the fact that 
the disk of the larger of the two progenitor spirals is not 
completely destroyed. 

(ii) The positive radial gradient in the velocity disper¬ 
sion of the young stars in E-|-A’s is not as pronounced in 


models with smaller m2. This is because there is less gaseous 
dissipation in a weaker starburst and thus the young stars 
with relatively cold kinematics are less centrally concen¬ 
trated in these models with smaller m2. Figures 19 and 20 
show the 2D distributions of velocity dispersion and the ra¬ 
dial profile of velocity dispersion, respectively, for old and 
young stars in the model Ul. It is clear from these figures 
that even the unequal-mass model shows a difference in kine¬ 
matics between the old and young stars. The variation in 
velocity dispersion a (km s“^) with radius, R (kpc) seen in 
these figures for R < 3 kpc can be described as: 

o- «-15.1 X 193.3 (16) 

for the old stars, and 

(7 « 0.2 X i?-b 90.9 (17) 

for the young stars. 

(iii) The mass ratio of young stars with EW(H5) > 2 A 
to all new (young) stars during the post-starburst phase, 
is significantly smaller for models with smaller m2 (< 0.1). 
This is because the secondary starburst is much weaker in 
such models, compared to the major or unequal-mass merger 
models. For example, the mass ratio of the strong-Hd young 
stars is only 13% in the model U2 with m2 = 0.1, and 58% 
for the fiducial model. This suggests that minor mergers 
with m2 less than 0.1 are highly unlikely to become E-|-A’s 
after their secondary starbursts, due to the very small frac¬ 
tion of A-type stars. The present models with m2 = 0.1 
show EW(H(5) of ~2 A and EW(H/I) of ^3.5 A in their post- 
starburst phases. Therefore, the observed disky E-|-A’s with 
strong Balmer absorption lines such as those observed in 
Zabludoff et al. (1996) are not likely to be formed from mi¬ 
nor merging. 

(iv) The formation efficiency of young GC’s depends 
on m2 in such a way that a larger number of GC’s are 
formed during merging in the models with larger m2. There¬ 
fore, E-bA’s with flattened (SO-like) morphologies formed 
via unequal-mass mergers show a smaller number of young 
GC’s than those with elliptical morphologies formed via ma¬ 
jor mergers. Therefore, the observed diversity in the number 
of young GC candidates in E-bA’s (Yang et al. 2004) could 
be due to this diversity of m2 in galaxy merging. Flattened 
E-bA’s are likely to show flattened spatial distributions of 
young GCs. 

3 . 2.4 Miscellaneous 

Finally, we note a number of other parameter dependences 
that our models have revealed, in particular ones involving 
/g, /b, and the IMF: 

(i) It is possible that two moderately strong starbursts, 
separated by a time interval of Gyr, can occur during 
major merging if the progenitor spirals have either no bulge 
(/b =0.0) or only a very small bulge (/b < 0.1). The first 
and stronger of these starbursts occurs as a result of the for¬ 
mation of a stellar bar, which in turn causes an inflow of gas 
into the center of the disk (e.g., in the bulgeless model M9). 
After this first starburst phase, the two galaxies are still 
separated yet at least one of them has a strong poststar- 
burst population. The second weaker starburst then occurs 
when the two disks finally merge to form an elliptical. Due 
to the rapid consumption of gas in the first starburst, the 
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Figure 20. The same as Fig. 3 but for the model T1 projected 
onto the x-y plane. 

final E+A elliptical will, after the second starburst, have a 
smaller fraction of A-type stars (39% in the M9 model, com¬ 
pared to 58% in the fiducial model) and thus have weaker 
H(5 absorption (EW(H5) = 5.3 A in the M9 model). In the 
interval between the two starbursts, such a merging system 
would appear as an interacting pair with an E-l-A spectrum, 
providing the star formation rate is insignificant. 

(ii) The Kinematical differences between the old and 
young stars are larger in the models with larger /g (i.e., 
larger gas mass fraction), essentially because a larger 
amount of random kinematical energy in the gas can be 
lost during merging, due to the more efficient gaseous dissi¬ 
pation in these models. The equivalent width of H5 during 
the E-l-A phases does not depend strongly on /g for a given 
IMF, as long as the merging galaxies are sufficiently gas-rich 

(/g > 0.1). 

(iii) The equivalent width of Hd depends strongly on the 
IMF of the starburst, in the sense that EW(H5) is larger for 
models with larger mu. For example, the simulated E-l-A in 
the Ml model with /g = 0.1 has EW(H5) = 7.4 A for Ml = 
IM 0 and EW(H(5)=1 .iA for Ml = 0.1 M©, whereas the 
E-l-A in the M7 model with /g = 0.5 has EW(H(5) = 6.9 A 
for Ml = IM0 and EW(H5) = 5.6 A for Ml = 0.1 M0. 
In contrast, EW(H/1) does not depend that strongly on the 
IMF. 


3.3 Tidal interaction 

In the tidal interaction model, a strong starburst will occur 
if the two disk galaxies are of comparable masses and they 
interact in a prograde manner i.e., if the intrinsic spin axis 
of the disk is parallel to the orbital axis. This is due to the 
formation of a stellar bar that drives gas into the central 
regions. Figure 20 shows the mass distributions of stars with 
different EW(H(5) in the tidal interaction model Tl, which 
develops a central stellar bar after the starburst. The mass 
distribution clearly depends on the EW(H(5) of the stars: 


H(5 
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Figure 21. The same as Fig. 5 but for the model Tl projected 
onto the x-y plane. The abscissa and the ordinate represent the 
x-axis and the i;-axis, respectively. 
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Figure 22. The same as Fig. 12 but for the model Tl projected 
onto the x-y plane. The abscissa and the ordinate represent the 
x-axis and the i;-axis, respectively. The grid points with no stars 
are shown in the darkest color (i.e., cr = 0kms“^). 

(i) the young stars (with EW(H(5) >6 A; bottom right-hand 
panel of Fig. 20) delineate a very compact bar in the center 
of the model, (ii) the old stars (with EW(H5) < OA; top 
left-hand panel of Fig. 20) delineate a global (~10kpc scale) 
bar with two spiral arms, and (iii) the intermediate age stars 
(with 0 < EW(H5) < 2A; top right-hand panel in Fig. 20) 
appear to have a ring-like morphology. 

Because of the elongated distribution of the Hd-strong 
stars, the 2D distributions of EW(H5) and EW(H/3) are also 
very flattened. Figure 21 clearly indicates that the radial 
gradients of EW(H(5) and EW(H/?) along the major axis 
of the flattened mass distribution of the young stars (with 
EW(H(5) > 6 A) is smaller than those perpendicular to the 
major axis, though the radial gradients are essentially neg¬ 
ative (i.e., large in the inner regions). The V — I and I — K 
colors also have 2D distributions with the same shapes and 
are also seen to have positive color gradients. Thus the sim¬ 
ulated E-l-A has positive color gradients and negative gra¬ 
dients in EW(H(5). 
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Figure 23. The same as Fig. 12 but for the tidal interaction 
model Tl, projected onto the x-z plane. The abscissa and the 
ordinate represent the x-axis and the 2 -axis, respectively. In order 
to more clearly see the kinematical difference between the young 
and old stars, the grid points with a < 20kms“^ (most of which 
are those where no young stars can be found) are shown in the 
darkest color. 
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Figure 24. The same as Fig. 13 but for the tidal interaction 
model Tl. The radial profile is derived for the 2D distribution 
shown in Fig. 23. 
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Figure 22 shows that the vertical velocity dispersion of 
the old stars is significantly larger than that of the young 
stars in the 2D distribution of velocity dispersion for the 
Tl model. This is because old stars are dynamically heated 
during tidal interaction whereas the gas, from which young 
stars form, dissipates its random kinematical energy during 
the interaction. The mean dispersion is 55.5kms“^ for the 
old stars and 28.8 kms“^ for the young stars. The radial de¬ 
pendence of the vertical velocity dispersion, a (km s“^), can 
be estimated from the 2D distributions shown in Figure 22 
and expressed as a function of R (radius from the center of 


the E-l-A) for R < 3 kpc: 


a ss —9.0 X R + 76.2 

(18) 

for the old stars, and 


a « -11.1 X 54.3 

(19) 


for the young stars. Thus the radial gradient in velocity dis¬ 
persion is steeper for the young stars than the old stars. 

Figures 23 and 24 reveal quite clear differences in the 2D 
velocity dispersion field (when viewed edge-on, i.e., in the or¬ 
bital plane of the interacting galaxy) between the young and 
old stars. Overall, the old stars show a higher velocity dis¬ 


persion than the young stars within the central 3 kpc of the 
simulated E-l-A, and a shallower radial gradient in the veloc¬ 
ity dispersion compared with the major merger models. The 
young stars show a monotonously increasing radial profile in 
velocity dispersion over the range -3 kpc< a: < 3 kpc along 
the major axis (i.e., a:-axis) of the mass distribution. This 
non-axisymmetric profile of velocity dispersion might be an 
indication that the young stars are in a non-equilibrium dy¬ 
namical state in the E-l-A phase. Such a non-axisymmetric 
profile in velocity dispersion cannot be seen in the 2D veloc¬ 
ity dispersion fields viewed from the x-axis and the a-axis, 
which implies that streaming of the gas (and young stars) 
along the stellar bar is partly responsible for the peculiar 
profile. 

Formation of stellar bars is essential for triggering 
strong starbursts and thus for producing young poststar- 
burst populations responsible for the E-l-A spectral sig¬ 
nature. Therefore, the models in which bar formation 
is strongly suppressed are less likely to become E-|-A’s. 
Amongst the tidal interaction models run for this study, T4 
(with the disk orbiting in a retrograde manner) and T5 (with 
the big bulge, fh — 1.0) are less likely to become E-|-A’s, due 
to the weak starbursts they experience during the tidal in¬ 
teraction. These results suggest that if E-|-A’s are formed 
from tidal interactions between gas-rich spirals, the E-|-A’s 
are more likely to be morphologically classified as barred 
disk galaxies, and more specifically as barred SOs due to 
the poor visibility of the spiral arms. Furthermore, a strong 
starburst associated with the formation of a stellar bar does 
not occur in the T6 model (where m2 is small, being only 
0 .1), because the tidal perturbation is not strong enough to 
form a bar. Only interaction models with a large m2 (> 0.3) 
show E-l-A spectra. This result provides a clue to the origin 
of the observed low luminosity disky E-|-A’s, because lower 
luminosity spirals are more likely to interact with galaxies 
more massive than themselves. Young GC’s are also formed 
in the interaction models, and have a ’disky’ spatial distribu¬ 
tion, quite different to that of the merger models described 
earlier. 


4 DISCUSSION 

4.1 Kinematical differences between old and 
young stars in E-|-A’s 

Norton et al. (2001) were the first to measure the spatially 
resolved kinematic properties of young and old stellar pop¬ 
ulations in E-|-A’s and they concluded that (i) both rapidly 
rotating and slowly rotating or non-rotating systems can be 
found in the E-l-A phase, (ii) most E-l-A galaxies (70%) show 
no evidence for rotation, (iii) the rotation is generally seen in 
both the old and young stellar populations in E-l-A galax¬ 
ies, (iv)the velocity dispersion ranges from ~30kms“^ to 
toto200kms“^, and (v) the young populations have on aver¬ 
age a higher velocity dispersion than the old populations. 
In addition, Norton et al’s measurements of the radial pro¬ 
files in velocity dispersion for their E-l-A galaxy sample (see 
their Fig. 4), leads us to conclude that (vi) about half of 
the galaxies show inwardly decreasing dispersion profiles for 
the young stellar populations, with some at least showing 
the young stars to have a smaller central velocity dispersion 
than the old stars (e.g., EA5, EA7, and EA15). 
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Conclusions (i)-(iii) above are consistent with our 
merger models in which E+A’s show a diversity in the de¬ 
gree of rotation of both the old and young stellar popula¬ 
tions. The fourth conclusion can also be understood in terms 
of differences in the masses of the merger progenitor spirals, 
because more massive E-|-A’s with a larger velocity disper¬ 
sion are formed in mergers between more massive spirals 
with larger rotational velocities. However, conclusion (v) ap¬ 
pears to be inconsistent with the present models in which 
the young stars in nearly all of the simulated E-|-A’s show a 
smaller central velocity dispersion than the old stars. This 
inconsistency might be due to the different methods that 
we and Norton et al. have employed to derive the velocity 
dispersion of the young stars. However, it could also be due 
to the fact that the interstellar gas in the simulations dis¬ 
sipates away much more random kinematic energy during 
galaxy merging than actually occurs. The dynamically cold 
young stellar populations in the simulations could well be 
due to ’over-dissipation’ in the gas. 

Gaseous dissipation, which is responsible for the rapid 
gas inflow that leads to a massive starburst during merg¬ 
ing, is essential for E-|-A formation as part of this process. 
Furthermore, the fact that our galaxy merger models are 
consistent with (vi) above also points to the need for effi¬ 
cient gaseous dissipation, since it is needed to produce dy¬ 
namically cold young stellar systems. Therefore, the incon¬ 
sistency between our models and conclusion (v) implies that 
the present numerical method for handling gaseous dissipa¬ 
tion and star formation does not allow a proper treatment 
of the interstellar gas dynamics in galaxy mergers. It there¬ 
fore remains unclear how the young stellar populations can 
have a larger velocity dispersion than the old stellar popu¬ 
lations, if they are formed from dissipative processes of star 
formation. 

Important factors other than the ’over-dissipation’ of 
gas could be responsible for the larger velocity dispersion 
of young stars in E-|-A’s. For example, if a massive black 
hole (MBH) and young A-type stars coexist in the central 
10 pc region of an E-l-A, the central velocity dispersion of 
the young stars will be significantly higher than that of the 
surrounding older stars, due to the deep potential well of 
the MBH. If this is the case, the observed diversity in the 
differences in central velocity dispersion between old and 
young stars in E-|-A’s could be ascribed to the existence or 
non-existence of MBHs in the center of E-|-A’s. Thus two key 
fines of investigation in future studies will be to investigate 
whether an improvement in our numerical modeling of the 
interstellar gas leads to more consistent agreement with (v) 
above, and whether central MBHs in the progenitor galaxies 
influence the radial profiles in velocity dispersion in E-|-A’s 
formed via interactions and merging. 

4.2 Origin of negative color gradients in E+A’s 

The color gradients in E-l-A galaxies can be considered to 
be one of the key physical properties which help us to deter¬ 
mine the most plausible physical mechanism(s) of E-l-A for¬ 
mation. This is because the spatial distribution of the most 
recent stars formed in an E-l-A depends strongly on how the 
star formation was truncated, and this is well traced by the 
color gradients (e.g.. Rose et al. 2001). Bartholomev et al. 
(2001) investigated the color gradients in 24 E-l-A galaxies, 


as well as those in 46 other galaxies with normal spectra. 
They found that, on average, the E-|-A’s had more positive 
color gradients than the normal galaxies. This is consistent 
with a merger/interaction origin for the E-l-A signature. 

However, the origin of E-l-A galaxies with negative color 
gradients - which are observed both in the field and in clus¬ 
ters (e.g., Bartholomew et al. 2001, Yang et al. 2004) - re¬ 
mains unclear. In particular, the central red colors in E-|-A’s 
with negative HS gradients (e.g., EA2 and EA4 in Yang et al. 
2004) are not easily explained within the interaction/merger 
scenario for E-l-A formation. We suggest here three possible 
explanations for the existence of E-l-A galaxies with negative 
color gradients: 

Firstly, the redder central colors maybe due to heav¬ 
ier dust extinction associated with the central poststarburst 
population. The question then is why do only a fraction of 
E-|-A’s have dusty gas surrounding their A-type stellar pop¬ 
ulations, and how can dusty gas still exist given that most 
of the residual gas should be blown away from the central 
regions of E-|-A’s by supernovae explosions associated with 
the secondary starburst? 

A second possible explanation is that the IMF during 
the starburst is truncated, in that no stars with masses 
smaller than 2 M© are formed. Hence the poststarburst pop¬ 
ulations would be dominated by red evolved stars (rather 
than by blue, main-sequence stars), which would give rise to 
the redder colors in the E-l-A phase (Chariot et al. 1993). 
The problem with this scenario is that it is unclear under 
what physical conditions the formation of stars with masses 
smaller than 2 Mq can be preferentially suppressed. 

Another possible explanation is that E-|-A’s with neg¬ 
ative color gradients are actually not poststarburst galax¬ 
ies but are starburst galaxies heavily obscured by dust. Al¬ 
though there is marginal evidence that some cluster galax¬ 
ies with E-l-A spectra show radio emission that could be 
powered by star formation activity (Smail et al. 1999), ra¬ 
dio continuum observations by Miller & Owen (2001) of the 
field E-|-A’s in Zabludoff et al’s (1996) sample yielded only 
2 detections, with radio luminosities consistent with only a 
moderate level of star formation (inconsistent with vigorous 
starbursts with 10 — 100 M©). Furthermore, ultra-luminous 
infrared galaxies, some of which are believed to be dusty 
starburst galaxies, do not have E-l-A spectra, but rather 
“e(a)” spectra where [OHJA3727 emission is seen in addi¬ 
tion to the strong Balmer line absorption (Poggianti & Wu 
2000 ). 

Unfortunately, the data currently available for E-l-A 
galaxies does not allow us to constrain the degree of dust 
extinction, the exponent of the IMF, nor the spatial distri¬ 
bution of dust. Therefore it is not possible at the moment 
to determine which of the above three explanations are the 
most likely. It should be stressed that we have assumed in 
this present study that most E-l-A galaxies are formed from 
galaxy interactions and merging. It might be necessary to 
relax this assumption to the extent that two or more phys¬ 
ical processes might be responsible for E-l-A formation. We 
will discuss this possibility in our forthcoming papers. 

So far we did not discuss color gradients of E-f As formed 
by merging between an old giant elliptical and a gas-rich 
spiral, because no theoretical studies (including the present 
one) have investigated spectrophotometric evolution of this 
type of galaxy mergers. Although it is not clear this ‘E-Sp’ 
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merging can create E+A’s, it is doubtlessly worthwhile for 
us to investigate this possibility in our future papers, given 
the fact that some Es show positive Balmer line gradients. 
E+A’s formed by E-Sp merging could show positive color 
gradients, if poststarburst components originally in the spi¬ 
ral progenitors are dispersed into the outer parts of the rem¬ 
nants (owing to the tidal destruction of the spirals) and thus 
located in the outer part of the remnants. 

4.3 Origin of kinematically decoupled cores in 
elliptical galaxies 

About one third of nearby luminous elliptical galaxies are 
observed to show peculiar core kinematics, such as counter- 
or oddly-rotating cores, and are thus considered to have 
rapidly spinning thick disks or torus-like components (see 
Bender 1996 for a review). Numerical simulations by Hern- 
quist & Barnes (1991) showed that ellipticals with counter¬ 
rotating cores can originate from dissipative major galaxy 
merging between two spirals. If E+A ellipticals are formed 
by dissipative major merging, as shown in the present study, 
the implication therefore is that a significant fraction of 
E+A’s must possess kinematical distinct cores (KDCs) that 
are very bright and have strong Hd absorption. Although 
Norton et al. (2001) have already investigated possible dif¬ 
ferences between the kinematics of the old and young stellar 
populations, they did not hnd signihcant kinematical differ¬ 
ences between the two populations. However, their results 
are based on long-slit spectroscopy of the stellar populations 
in only the central few kiloparsecs of E+A galaxies, and it 
remains unclear whether the kinematics of this region are 
different to those in the more outer regions. 

The present simulations have demonstrated that an 
E+A galaxy formed in a major merger will have a very flat¬ 
tened 2D distribution in H5, resulting from the central disky 
core composed mostly of A-type stars. They have also shown 
that the 2D kinematical distributions of the young stars are 
significantly different from those of the old stars, which re¬ 
flects the fact that the young stellar population in the core is 
more strongly supported by rotation. These results suggest 
that E+A galaxies might well be the best ’signposts’ of KDC 
formation in dissipative major merging. Future spatially re¬ 
solved, integral held unit spectroscopy of E+A galaxies with 
8-10 m class telescopes should be capable of mapping their 
2D kinematical and spectral properties to the level required 
to address the origin of KDCs in ellipticals. If such observa¬ 
tions also confirm the diversity in the 2D Hd and kinematical 
distributions seen in our simulations, they will also conhrm 
the dissipative major merger scenario for E+A formation. 


4.4 Formation of disky E+A’s 

Recent morphological studies of E+A’s have revealed that 
they are not all spheroidal systems, with a certain fraction 
having disks, and some of these being disk-dominated (Blake 
et al. 2004; Tran et al. 2004). Moreover, there seems to be 
environmental differences in that more ’disky’ E+A galax¬ 
ies are seen in rich clusters (up to ~20% of the E+A pop¬ 
ulation; Tran et al.) than in the lower-density group and 
held environment (only a ^5% fraction; Blake et ah). These 
trends might suggest that different physical mechanisms are 


responsible for E+A formation in clusters and in groups/the 
held. Ram-pressure effects in clusters of galaxies are sug¬ 
gested to enhance star formation activity in gas-rich spi¬ 
rals without destroying their disk components (e.g., Bekki 
& Couch 2003, Milvang-Jensen et al. 2003; Bamford et al. 
2004) and thus can be considered to contribute to the large 
fraction of disky E+A’s in clusters. Since galaxy interac¬ 
tions and merging are thought to only take place in the 
galaxy group and held environments, the presence of disky 
E+A’s raises the question as to what kinds of galaxy inter¬ 
actions and merging lead to the formation of E+A’s with 
this morphology. 

Our numerical simulations have demonstrated that: 
(i) rotationally supported and battened E+A’s with an SO 
morphology, can be formed in some unequal-mass mergers, 
and (ii) E+A disk galaxies with thick disks can be formed 
in galaxy interactions with special initial parameters (e.g., 
prograde encounters). The latter of these is always asso¬ 
ciated with bar formation and double spiral arms (rather 
than multiple-arm structures). These results suggest that 
unequal-mass mergers are only responsible for the formation 
of E+A’s with an SO morphology, and that tidal interactions 
can be responsible for the formation of E+A’s with barred 
morphologies or those with double-arm structures. They also 
suggest that held E+A’s without any thick disk component 
can only be formed by galaxy merging and interaction. The 
present numerical simulations have also shown that the ver¬ 
tical velocity dispersion of old stars in disky E+A’s can be 
rather high (up to 60 km s“^) owing to the dynamically hot 
thick disk. Therefore high-resolution imaging of the vertical 
structure in disky E+A’s and spectroscopic determination of 
the velocity dispersion of their old stars can provide valuable 
information on the presence of dynamically hot thick disks 
and thus assess the viability of the disky E+A formation 
scenario described in this paper. 


5 CONCLUSIONS 

We have investigated the dynamical and spectrophotomet- 
ric properties of E+A galaxies formed via galaxy interac¬ 
tions and merging, using gas dynamical simulations com¬ 
bined with stellar population synthesis codes. Our principle 
results can be summarised as follows: 

(1) E+A ellipticals formed by dissipative major galaxy 
merging show positive radial gradients in color (i.e., bluer 
colors at their center) and negative radial gradients in 
Balmer absorption line strength, due to the larger fraction 
of A-type stars in their inner regions. These color and line 
index gradients become shallower as time passes by, because 
of the aging of the poststarburst stellar population. These 
numerical models, however, cannot explain the E+A galax¬ 
ies that are observed to have negative radial gradients in 
both color and Balmer absorption line strength. 

(2) The dynamical and spectroscopic properties of E+A 
ellipticals formed by dissipative major galaxy merging have 
2D distributions which can be remarkably different, depend¬ 
ing on the orbital parameters of the merger. Furthermore, 
such differences also exist between the old and young stars 
in these galaxies. For example, E+A ellipticals with kine¬ 
matically decoupled cores (KDC’s) have a very flattened HJ 
absorption distribution, with differences in rotation and ve- 
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locity dispersion between the old and young stars. Future 
spatially resolved, integral field unit spectroscopy with 8- 
10m class telescopes should be able to bring the past KDC 
formation sites in young E+A ellipticals into relief and thus 
provide an evidence that KDCs can be formed from dissipa¬ 
tive major merging. 

(3) The 2D distributions of colors and Balmer absorp¬ 
tion lines in E+A ellipticals show a larger internal dispersion 
compared with those of ’passive’ ellipticals with weak H(5 
(~0A). The internal color dispersion of ellipticals formed by 
major merging becomes smaller as the E+A phase passes, 
and therefore there will be an anti-correlation between the 
size of this dispersion and H(5 strength in the post-E+A 
phase. These results imply that the large color dispersion in 
the 2D photometric properties of E+A ellipticals provides 
further evidence for the transformation from gas-rich late- 
type spirals into passive ellipticals. 

(4) Unequal-mass galaxy merging (m2 ~ 0.3) can trig¬ 
ger central starbursts and thus can transform two gas-rich 
spirals into very flattened ellipticals or SOs with an E+A 
spectral signature . These E+A’s are flattened by rotation, 
with a large V/a (> 1), and show positive radial gradients 
in color and negative radial gradients in Balmer line ab¬ 
sorption. These flattened E+A’s could be the missing link 
between gas-rich late-type spirals and passive SOs. 

(5) The mass fraction of young stars with H(5 > 2 A 
in a merger remnant strongly depends on the mass ratio 
(m2) of the two progenitor spirals in the sense that the frac¬ 
tion is larger in merger remnants with larger m2. This is 
essentially because secondary starbursts during merging are 
stronger and thus produce a larger fraction of A-type stars. 
Therefore, it is highly unlikely that minor mergers with 
m2 < 0.1 are responsible for E+A formation with strong 
H(5 (EW> 6 A). 

(6) Strong tidal interaction between galaxies can also be 
responsible for E+A formation, although special orbital con¬ 
figurations (i.e., prograde encounters, larger m2, and smaller 
bulges) are required for the tidal encounters. The E+A’s 
formed from tidal interactions are likely to have an SBO mor¬ 
phology, with thick disks, positive color gradients and neg¬ 
ative H(5 gradients. Our results thus suggest that some frac¬ 
tion of disky E+A’s with SBO morphology may be formed 
via strong tidal interactions. 

(7) E+A’s formed by galaxy merging and interactions 
can have young, bright, metal-rich, and H(5-strong GC’s, 
mostly in the central few kiloparsecs, because the number of 
cloud-cloud collisions with high relative velocities and small 
impact parameters are dramatically enhanced during this 
process. Disky E+A’s are more likely to show a disky distri¬ 
bution of young GC’s and rotational kinematics in the GC 
population. Eurthermore, the total number of young GC’s 
is more likely to be larger in E+A ellipticals than in disky 
E+A’s for a given luminosity, because their formation effi¬ 
ciency is higher in major mergers than in minor mergers and 
tidal interactions. 

(8) The present models of E+A formation have difficul¬ 
ties in explaining: (i) E+A’s showing both negative radial 
gradients in color and EW(H(5), and (ii) E+A’s where the 
central velocity dispersion of the young stars is larger than 
for the old stars. These difhculties suggest that more so¬ 
phisticated numerical methods for star formation and inter¬ 
stellar dynamics are required to reproduce self-consistently 


the diversity in the radial gradients of dynamical and spec- 
trophotometric properties of the central poststarburst pop¬ 
ulations in E+A’s. 
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